HD-THEP-94-35 



Kosterlitz-Thouless phase transition 
in the two dimensional linear cr-model 

M. Grater and C. Wetterich 

Institut fiir Theoretische Physik, Universitdt Heidelberg, 
Philosophenweg 16, 69120 Heidelberg 



Abstract: 

We investigate the 0{N) symmetric linear cr-model in two dimensions by means of an exact 
nonperturbative evolution equation. The perturbative infrared divergences are absent in this 
formulation. We use a simple approximative solution of the flow equation which corresponds 
to a derivative expansion for the effective action. For N = 2 this gives a good picture of the 
Kosterlitz-Thouless phase transition. 



In dimensions smaller than four the A'^-component linear cj-model is not directly accessible 
to perturbation theory. The obstacle arises from infrared divergences due to massless modes 
- either the Goldstone modes in the regime with spontaneous symmetry breaking for N > 2 
or the massless mode at the phase transition. One usually ressorts to an extrapolation of the 
4 — e expansion around the four dimensional model, or to the 2 + e expansion around the 
corresponding two dimensional nonlinear (j-model. We present here an approach working in 
arbitrary dimension d which can cope with the infrared problems. It is based on the concept of 
the average action Q which is the continuum analogue of the block spin action on the lattice 

The computation of includes all fluctuations with momenta q"^ larger than the infrared 
cutoff k'^. For A; ^ one recovers the effective action, i.e. the generating functional for the one 
particle irreducible Green functions or the free energy. The dependence of on the scale k is 
determined by an exact evolution equation 0]: 



|r. = irr{(rr + . (i) 



dt 

(In a Fourier basis the trace reads Tr = Y,a=i I i^J^d t = ln(/c/A).) In distinction to earlier 
versions of exact renormalization group equations the flow equation (||) is characterized 
by its close resemblance to perturbation theory: The main differences as compared to a one 
loop expression are the appearance of the exact inverse propagator F^ (the second functional 
variation of F^ with respect to the A^-component scalar fields and the infrared cutoff R^- 
With 

2n Zkq^ exp(-gVfc^) 

Rk[q ) = 7 — TrPK 

1 — exp(— g^//c^j 

(where is an appropriate wave function renormalization constant specified below) we observe 
that the momentum integral in (0) is both infrared and ultraviolet finite. The short distance 
physics has to be specified by the "initial value" Fa (the "classical action") at some high mo- 
mentum scale A. A solution of the flow equation (|^) allows to extrapolate from the classical 
action Fa to the effective action Fq and therefore constitutes a solution of the model. 

Our nonperturbative evolution equation remains a complicated nonlinear functional differen- 
tial equation which cannot be solved exactly. One needs approximate nonperturbative solutions 
which correspond to a truncation of the general form of F^ and therefore reduce the problem 
to a manageable number of degrees of freedom. We consider here the 0{N) symmetric linear 
o"- model with a truncation corresponding to the lowest terms in a derivative expansion of F^, 
i.e. 

Ffc = 1 (fx \uk{p) + ^Zkd^'<Padf,r] , (3) 

where p = ^4>a(f>°'- The truncation (^ neglects all terms with more than two derivatives. The 
most general two derivative expression would contain for > 1 a second term Y^d^pd^p and 
Zk and Yk could depend on p. For simplicity we consider first only a constant Z^ and Y^ = 0. 

With the Ansatz (P) we obtain from (jl]) an evolution 

equation for the effective average potential: 

-C/.W=../ ...*^»,(.)(_ + _) (4) 



where x = q with 



Mo = xil + n{x)) + Z^'uUp) (5) 
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Ml = x{l + rk{x)) + Z-^U',ip) + 2Z-^pU'^{p) 

Here primes denote derivatives with respect to p and we observe the appearance of p depen- 
dent mass terms ~ C/^ in the "Goldstone" and ~ + 2pC/^' in the "radial" directions. The 
dimensionless functions Sk and depend on the ratio x/k'^ and are given by 



Rkjx) 
Z^x 



The second term in Sk{x) is proportional to the anomalous dimension 

d 

'? = -^ln2'fc (7) 

and arises from the factor Zk in @. We will neglect it in the following. 
It is convenient to work with dimensionless and renormalized quantities: 

p = k^-''Zkp ; Uk{p) = k-''Uk{pip)) y = ^- (8) 

To evaluate the equation for the potential we make a further approximation and expand around 
the minimum of u/. for non zero p = k up to the quadratic order in p: 

Uk{p) = Uk{K) + ^X{p - Kf . (9) 

The condition du/dp\p=K = holds independent of t and gives us the evolution equation for the 
location of the minimum of the potential parametrized by k,. A similar evolution equation can 
be derived ||5| for the symmetric regime where k = and an appropriate variable is du/dp\p=o. 
The flow equations for k and A follow from @) by variable substitution and appropriate differ- 
entiations with respect to p: 

P. = ^ = -{d-2 + ri)K + 2vd{N-l)lf{0)+6vdli{2XK) 

Px = = id-4 + 2r])X + 2vdiN -l)X'^4{0) + 18vdX'^4{2XK). (10) 

We notice that the "threshold functions" 

liiu) = -J^ dyy-2s{yMl + r{y))+u]-(-^'^ (11) 

vanish for u; — > oo and describe appropriately the decoupling of a heavy radial mode for 2Xk 3> 1. 
The integrals are normalized such that (0) = 1, /^(O) = 1- -^l- ( |lO|) describes how the average 
potential Uk changes its shape, with "initial condition" given at some short distance scale k = A 
by k(A) and A(A). The solution for k = determines the IPI two and four point functions at 
zero momentum. In particular, the theory is in the symmetric phase if «;(0) = - this happens 
if K reaches zero for some nonvanishing kg > 0. On the other hand, the phase with spontaneous 
symmetry breaking corresponds to po{0) > where po{k) = k'^^^ Z^^^ K(k). A second order phase 
transition is characterized by a scaling solution corresponding to fixpoints for k. and A. For small 
deviations from the fixpoint there is typically one infrared unstable direction which is related to 
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the relevant mass parameter. The phase transition can be studied as a function of k(A) with a 
critical value ac(A) = Kc- The difference k(A) — Kc can be assumed to be proportional to Tc — T, 
with Tc the critical temperature. This allows to define and compute critical exponents in a 
standard way. We should mention a particular possibility for d = 2, namely that k{0) remains 
strictly positive whereas po{0) vanishes due to limfc^o c«- This is a somewhat special form 

of spontaneous symmetry breaking, where the renormalized expectation value, which determines 
the renormalized mass, is different from zero whereas the expectation value of the unrenormalized 
field vanishes. We will see that this scenario is indeed realized for d = 2, N = 2. The phase with 
this special form of spontaneous symmetry breaking exhibits a massive radial and a massless 
Goldstone boson - and remains nevertheless consistent with the Mermin- Wagner theorem Q that 
the expectation value of the (unrenormalized) field (pa niust vanish for N > 2. The Kosterlitz- 
Thouless phase transition ||^] describes the transition from this phase to the standard symmetric 
phase of the linear cj-model, i.e. the phase where k{0) = with a spectrum of two degenerate 
massive modes. 

In order to solve the flow equation ( [l0| ) we further need the anomalous dimension t] which 
is related to the scale dependence of the wave function renormalization. We identify the wave 
function renormalization in the infrared cutoff (||) with the coefficient multiplying the kinetic 
term in the average action (Q), evaluated for zero momentum and p = k. The flow of Zk can then 
be computed from the momentum dependence of -§1^^^^ Q and one obtains with the truncation 
® 

r) = ^X\ml2(.0,2XK). (12) 
The integral m2 2 defines again an appropriate threshold function 

.(13) 



"^2,2(0,w) = -^y^ dyy'^'\A\^ + r{y)+yQ^r{y)^ [y(l+r(y))]-2[y(l+r(y)) + 



Here ^ stands symbolically for a derivative which acts only on the infrared cutoff Rk in r, with 

^r[y) = s{y), and we omit again the second term ~ 77 in s We have now a system of 
three coupled non linear differential equations for k, A and rj which can be solved numerically 
for arbitrary and d. 

We next specialize to the two dimensional linear c-model {d = 2). We observe that the 
flow equations can be solved analytically in the limiting case of a large mass oo = 2Xk of the 
radial mode. The threshold functions vanish with powers of u>~^ and for > 1 the leading 
contributions to the /3-functions are those from the Goldstone modes. Therefore this limit is 
called the Goldstone regime. In this approximation the /3-functions can be expanded in powers of 
uj^^. In particular, the leading order of the anomalous dimension can be extracted immediately 
from (|12|): 

, = ji: + 0(.-). (14) 



Inserting this result in (llQ) we have 



= + 0{^~') (15) 



and the leading order of /3a is 



/?. = -2A + i^^-in^^A^ + 0(O. (16) 
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Eq. dH) has a fixpoint solution A* = (^j^ ^1)1^2 ^ 18.13/(iV - 1). 

For N > 2 there exists a simple relation between the linear and the nonlinear cr-model: The 
effective coupling between the Goldstone bosons of the nonabelian nonlinear cr-model can be 
extracted directly from (^) and reads in an appropriate normalization 

The lowest order contribution to P^. (0) coincides with the one loop expression for the running 
of as computed in the nonlinear a-model. We emphasize in this context the importance of the 
anomalous dimension rj which changes the factor (A^ — 1) appearing in ( p!0| ) into the appropriate 
factor (A^ — 2) in (|l5|). In correspondence with the universality of the two loop /3- function for 
in the nonlinar cj-model we expect the next to leading term ~ in (15) to be also 



proportional to {N — 2). In order to verify this one has to go beyond the truncation (y) and 
systematically keep all terms contributing in the appropriate order of k~^. (This calculation is 
similar to the extraction of the two loop /3-function of the linear cr-model in four dimensions 
by means of an "improved one loop calculation" using the flow equation (|l|) |^].) We have 
calculated the expansion of (3^ up to the order 0{k^^) for the most general two derivative 
action, i.e. neglecting only the momentum dependence of and Y^. The result agrees with the 
two loop term of the nonlinear cr-model within a few per cent, and the discrepancy should be 
attributed to the neglected momentum dependence of the wave function renormalization. The 
issue of the contribution to (3^ in order is less clear: Of course, the direct contribution of 
the Goldstone bosons (combined with their contribution to rj) should always vanish for N = 2 
since no nonabelian coupling exists in this case. The radial mode however, could generate a 
contribution which is not proportional to (A'^ — 2). This contribution is possibly nonanalytic in 
and would correspond to a nonperturbative contribution in the language of the nonlinear 
cr-model. 

Let us now turn to the two dimensional abelian model (d = 2, N = 2) for which we want to 
describe the Kosterlitz-Thouless phase transition. In the limit of vanishing for large enough k 
the location of the minimum of Uk{p) (^) is independent of the scale k. Therefore the parameter 
K, or, alternatively, the temperature difference Tc—T, can be viewed as a free parameter. If we go 



beyond the lowest order estimate (16) the fixpoint for A remains, but A* becomes dependent on 
K. This implies that the system has a line of fixpoints which is parametrized by k as suggested by 
results obtained from calculations with the nonlinear cr-model ||lO[]. In particular, the anomalous 



dimension r] depends on the temperature Tc — T (14). Even if this picture is not fully accurate 
for nonvanishing /3k, it is a very good approximation for large k: The possible running of k is 
extremely slow, especially if /3k vanishes in order k~^. We associate the low temperature or large 
K phase with the phase of vortex condensation in the nonlinear cr-model. The correlation lenght 
is always infinite due to the Goldstone boson. Since rj > we expect the inverse propagator 
of this Goldstone degree of freedom ~ ((7^)^ thus avoiding Coleman's no go theorem [|l^] 
for free massless particles in two dimensions. On the other hand, for small values of Xk the 



threshold functions (11)([T^) can be expanded in powers of Xk. The anomalous dimension is 
small and k is driven to zero for A;^ > 0. This corresponds to the symmetric phase of the linear 
cr-model with a massive complex scalar field. We associate this high temperature phase with 
the phase of vortex disorder in the picture of the nonlinear cr-model. The transition between 
the behaviour for large and small k is described by the Kosterlitz-Thouless transition. In the 
language of the linear cr-model it is the transition from a special type of spontaneous symmetry 
breaking to symmetry restoration. 

Finally we give a summary of the results obtained from the numerical integration of the 
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evolution equations (10) and ( [T^ ) for the special case d = 2, N = 2 ^ We use a Runge-Kutta 
method starting at t = with arbitary initial values for k and A and solve the flow equations 
for large negative values of t. For the integrals defined by ([Tl|) and ( p!3|) we use numerical fits. 
Results are shown in figs. 1-4 where we plot typical trajectories. The distance between points 
corresponds to equal steps in t such that very dense points or lines indicate the very slow running 
in the vicinity of fixpoints. 

The understanding of the trajectories needs a few comments: The work of Kosterlitz and 
Thouless 1^ suggests that the correlation lenght is divergent for all temperatures below a critical 
temperature Tc and that the critical exponent rj depends on temperature. The consequence for 
our model is that above a critical value for /t all (3- functions should vanish for a line of fixpoints 
parametrized by k,. From the results in the Goldstone regime and from earlier calculations [|l^] 
we conclude that (3^. should vanish faster than k^^ for large k. Our truncation (^), however, 
yields a function which vanishes only like k^^. The consequence is that even if the system 
reaches the supposed line of fixpoints the parameter k decreases very slowly until the transition 
to the symmetric regime is reached. The anomalous dimension first grows with decreasing k 
(14) until the critical value is reached. Then the system runs into the symmetric regime and r] 
vanishes. So we expect that rj reaches a maximum near the phase transition. We use this as a 
criterion for the critical value Kc- In summary, the truncation (|3|) smoothens the phase transition 
and this prevents a very accurate determination of the critical value Kc and the corresponding 
anomalous dimension rjc- From the numerical point of view the absence of a true phase transition 
in the truncation (P) makes life easier: One particular trajectory can show both the features of 
the low and the high temperature phase since it crosses from one to the other. 

The numerical results fullfil our expectations. Fig.|^ shows the evolution of the anomalous 
dimension with decreasing t for several different initial values k(A),A(A). The maximum is 
reached with rjc = 0.24 which has to be compared with the result of Kosterlitz and Thouless 
rjc = 0.25 IQ]. The approximate "line of fixpoints" for k, > Kc {v < Vc) is demonstrated by the 
selfsimilarity of the curves for large —t. Trajectories with different initial conditions hit the 
line of fixpoints at different k. Subsequently they follow the line of fixpoints. Except for the 
value of K{k) all "memory" of the initial conditions is lost for t < —3. Another manifestation 
of the line of fixpoints in the (k. A) plane is demonstrated in fig.^. After some fast "initial 
running" (dotted parts of the trajectories) all trajectories with large enough k{A) follow this 
line independent of the initial A(A). We emphasize that the nonlinear cj-model corresponds to 
A(A) oo. Our investigation shows that the linear u-model is in the same universality class, 
even for very small A(A). In figj^ we plot r/(K). Along the line of fixpoints we find perfect 
agreement with the analytical estimate (|l^) for large k. For k = 0.3 the deviation from the 
lowest order result is 34% in the present truncation. Finally we show in fig.§ the value of P\ for 
different trajectories. The dense parts of the "ingoing curves" show the fixpoint behaviour at 
A*(ki) where ki denotes the value of n where the line of fixpoints is hit. (For small k(A) there is 
a substantial difference between Kiand k(A) which depends also on A(A). This can be seen from 
the curves with k(A) = 1.) After hitting the line of fixpoints the trajectories stay for a large 
t-interval at I3x very close to zero. Subsequently, the "outgoing curve" indicates the transition 
to the symmetric phase. 

In conclusion, both the analytical and the numerical investigations demonstrate all impor- 
tant characteristics of the Kosterlitz-Thouless phase transition for the linear cr-model. This 
belongs to the same universality class as the nonlinear cr-model and we have demonstrated a 
close correspondence between the linear and the nonlinear cr-model with abelian symmetry. In 

^This work has also been done for d — 2, N = 1. There we find a fixpoint which corresponds to the second 
order phase transition in the Ising modeL 
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particular, the phase with vortex disorder in the nonhnear cj-model corresponds simply to the 
symmetric phase of the linear cr-model. We emphasize that we have never needed the explicit 
investigation of vortex configurations. The exact nonperturbative flow equation includes au- 
tomatically all configurations. Its ability to cope with the infrared problems of perturbation 
theory is confirmed by the present work. 

Despite the simple and clear qualitative picture arising from the truncation (^) this letter only 
constitutes a first step for a quantitative investigation. It is not excluded that the coincidence 
of our critical r]c ~ 0.24 with 1/4 is somewhat accidental. In order to answer this question one 
needs to go beyound the truncation (^). In view of the relatively large value of r]c we expect 
in particular that the momentum dependence of the wave function renormalization (or the 
deviation of the inverse propagator from g^) could play an important role at the phase transition. 
This effect should be included in a more detailed quantitative investigation. 
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Figure Captions 

• fig 1: The scale dependent anomalous dimension rj 

• fig 2: Critical line in the k — \ plane 

• fig 3: The k dependence of the critical exponent rj 

• fig 4: The /3-function for A 
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